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چکیده 


در بسیاری موارد بعد از وقوع سیل به اطلاعاتی درباره مشخصات سیلاب در بالادست یک موقعیت مشخص نیاز است و ممکن است رودخانه در آن 


نقطه فاقد ایستگاه آب‌سنجی باشد یا داده‌برداری در زمان سیلاب انجام نشده باشد. در روندیابی معکوس محاسبات از مقطع پایین‌دست به صورت گام به 
گام به سمت مقطع بالادست رودخانه صورت می‌گیرد و هیدروگراف ورودی بر اساس مشخصات هیدرولیکی رودخانه و هیدروگراف پایین‌دست تولید 
در خشک رودها بسیار مهم است. با جفت کردن رابطه نفوذ گرین- آمپت در معادله پیوستگی جریان انجام شده است. باوجود ماهیت ناپایدار مدل‌های 


روندیابی معکوس: به ویژه زمانی که عمق یا دبی اولیه جریان صفر است؛ با انتخاب پارامترهای وزنی مناسب ورودی به مدل» این فرایند با موفقیت انحام 
شد. قابلیت کاربرد مدل عددی توسعه یافته با استفاده از دو سری داده اندازه‌گیری شده و واقعی بررسی شد و نتایج نشان داد مدل عددی پیشنهادی. دبی 
اوج هیدروگراف ورودی يا بالادست را تا دقت ٩٩‏ درصد و زمان اوج را تا دقت ٩۷‏ درصد بازیابی نموده است. این ننایج قابلیت استفاده از مدل‌های 


روندیابی هیدرولیکی معکوس سیلاب در خشک‌رودها را تایید می‌نماید. 


واژه‌های کلیدی: روندیابی معکوس, ناپایداری مدل عددی» محاسبات گام به گام مقدار نفوذ معادله گرین- آمپت 


مقدمه 


رود و بستر آن» یعنی رودخانه نقش مهمی در زندگی بشر دارند. 
سیلاب‌های مهیب در رودخانه‌هاء تهدیدی برای تاسیسات» مناطق 
شهری. صنعتی و کشاورزی مجاور آن محسوب می‌شود و سالانه 
باعث ایجاد تلفات جانی و خسارات مالی فراوان می‌شود. 

برای جلوگیری و یا کاهش آثار مخرب سیل, شناخت ویژگی‌های 
سیل شامل ارتفاع آب و دبی متناظر با آن به‌ویژه دبی اوج سیلاب با 
استفاده از روش‌های روندیابی حائز اهمیت آنتست: روندیابی سیلاب از 


۱ و ۲- به‌ترتیب کاندیدای دکتری سازه‌های آبی و دانشیار گروه مهندسی آبیاری و 
آبادانی» دانشگاه تهران 
(#- نویسنده مسئول: 
۳- دکتری در مهندسی عمران» رئیس گروه تحقیقات کاربردی. شرکت آب 
منطقه‌ای خراسان رضوی 

(۵ ۷ 6 


6۵006 1.2.100912۷1ظ :1تمحصط) 


دو دیدگاه هیدرولوژیکی و هیدرولیکی قابل مطالعه است. روش‌های 
هیدرولیکی بر اصول تلوری جریان‌های غیردائمی در آبراهه روباز و 
ص ادا کال مطه ان سانشان حاشان اسف 
روش‌ها دقیق ولی روابط پیچیده و استفاده از آن‌ها مشکل است. در 
روش‌های هیدرولوژیکی, اصل پیوستگی جریان و رابطه بین دبی و 
ذخیره آب مورد استفاده قرار می‌گیرند. این مساله سبب می‌شود 
استفاده از این روش‌ها نسبت به روش‌های هیدرولیکی» حل مسائل 
روندیابی را به میزان قابل ملاحظه‌ای ساده نمایده ولی جواب‌های به 
دست آمده دارای دقت کمتری نسبت به روش‌های هیدرولیکی هستند 
و۱۸ 

تعیین هیدروگراف مقطع پایین‌دست با استفاده از هیدروگراف 
ورودی یا بالادست. روندیابی (یا روندیابی مستقیم سیلاب) نامیده 
می‌شود. در روندیابی معکوس» محاسبات از مقطع پایین‌دست به 
صورت گام به گام به سمت بالادست مسیر جریان آب انجام می‌شود 
و هیدروگراف ورودی (در بالادست یک بازه) با توجه به مشخصات 


۶ نشریه آب و خاک. جلد ۰۳۵ شماره ۶ بهمن - اسفند ۱۴۰۰ 


هیدرولیکی رودخانه و هیدروگراف در پایین‌دست تعیین می‌شود (۱۳). 
مشخصات سیلاب رخ داده در نقطه بالادست یک ایستگاه 
هیدرومتری نیاز است و ممکن است رودخانه در آن نقطه فاقد ایستگاه 
هیدرومتری باشد و یا داده‌برداری در زمان سیل اتفاق نیفتاده باشد. 
هدف از روندیابی معکوس بعد از مشخص کردن دبی یا تراز سطح آب 
سیلاب در بالادست» برآورد ابعاد مناسب سازه‌های مرتبط با م‌هنداسی 
است. از دیگر موارد کاربرد روندیابی معکوس می‌توان به واسنجی 
ایستگاه‌های هیدرومتری احداث شده در بالادست اشاره کرد (۲). 
روندیابی معکوس در برنامه‌ریزی بهره‌برداری از دریچه‌های تخلیه 
مخازن و سدهای مخزنی در مواقع سیلاب. برای جلوگیری از خسارت 
به مناطق پایین‌دست نیز قابل استفاده است (۱۲). 

نتایج بررسی مطالعات انجام شده در زمینه روندیابی معکوس 
سیلاب بیانگر این موضوع است که روش‌های هیدرولیکی که به حل 
معادله سنت- ونانت می‌پردازند و در صدر آن‌ها روش موج دینامیکی» 
8 ۰ ۱۱ و9 ۲ 

همچنین با توجه به دشواری کاربرد روش‌های هیدرولیکی در 
روندیابی معکوس سیلاب» محققان مختلف از روش‌های ساده‌تر 
استفاده کرده‌اند که هم به لحاظ زمان و هم به لحاظ هزینه به 
ماسکینگام و الگوریتم ژننیک از جمله روش‌های به کار رفته در 

داس (۱۱) از مدل ماسکینگام برای روندیابی معکوس سیلاب 
استفاده کرد. در این تحقیق» روشی برای تعیین پارامترهای معادلات 
ماسکینگام ارائه شد که برای حل غیرخطی معادلات روندیابی 
معکوس کاربرد داشت. نتایج نشان داد که برای روندیابی معکوس 
جریان, انجام واسنجی پارامترهای مدل الزامی است. 

کوسیس و همکاران (۱۹) به بررسی روندیابی معکوس سیلاب با 
هدف رون دیابی جریان ورودی از روی هیدروگراف خروجی یا 
پایین‌دست با استفاده از روش ماسکینگام پرداختند و نتیجه را با حل 
معکوس معادلات سنت- ونانت مقایسه کردند. نتایج آن‌ها دقت قابل 
قبول روش ماسکینگام را نسبت به روندیابی هیدرولیکی معکوس 

بادفر و همکاران (۲) عملکرد مدل ماسکینگام خطی و غیرخطی 
را در روندیابی معکوس هیدرولوژیکی سیلاب از طریق حل تابع ذخیره 
به روش اویلر و رانج کوتای مرتبه چهار مورد بررسی قرار دادند. 
(565) بهینه گردید. نتایج تحقیق آن‌ها در بازیابی هیدروگراف مقطع 
بالادست جریان دقتی بالغ بر ۸۲ درصد را در شرایط استفاده از 


داده‌های واقعی نشان داد که بیان گر قابلیت کاربرد مدل ارائه شده در 
شرایط واقعی جریان‌های رودخانه‌ای می‌باشد. 

در هنگام وقوع سیل در رودخانه‌هاء در اثر تبخیر و نفوذ از بستر و 
کناره‌هاء حجم سیلاب به مرور کاهش می‌یابد. این کاهش حجم 
نلاب که لفات اتقال تامیتهم شوخ در اطی شک قایل ترجه 
بوده و در برخی از مطالعات مقدار آن تا۲۰۱ درصد حجم سیلاب 
ورودی به بزه گزارش شده است (۶) از این‌رو مدل‌های توسعه یافته 
برای روندیابی سیلاب در رودخانه‌های فصلی و رودخانه‌های جاری در 
مناطق خشک و نیمه خشک باید بتوانند برآوردی مناسب از تلفات 
خصوص رواناب سطحی و تغذیه آب زیرزمینی به‌دست می‌دهد. مورین 
9 همکاران (۲۲) با انحام روندیابی سیلاب 9 با استفاده از داده‌های 
مشاهداتی سیلاب های رودخانه کوسیب در کشور نامیبیا که در منطقه 
بسیار خشک آفریقا واقع شده است بیان داشتند سیلاب‌های با حجم 
جریان متوسط تا بزرگ نقش مهمی در تغذیه سفره‌های آب زیرزمینی 
داشته که مقدار نشت در بسترهای با نفوذپذیری بالا و مجاری آبی 
طولانی و عریض محسوس‌تر است. شدت نفوذ در برخی وقایع 
سیلاب به حدی بود که در عمل هیچ‌گونه جریانی به مقاطع 
پایین‌دست رودخانه نمی‌رسید. 

اگر چه اهمیت تلفات انتقال در مناطق خشک از سالیان دور مورد 
توجه بوده است اما اطلاعات چندانی از فرآیند حاکم بر شکل‌گیری 
این تلفات در دست نیست. به علاوه. وجود این رودخانه‌ها در مناطق 
خشک و محروم که اغلب فاقد ایستگاه های اندازه گیری می‌باشند. 
برآورد تلفات انتقال در این مناطق را دشوار می‌نماید (۴). 

در کشور ما ایران» حدود ۳۹ درصد یعنی ۲۶,۰۰۰ کیلومتر از طول 
ریدشانه‌هاس اضای میم را وهای فصا شک کین 
داده‌اند ولی متاسفانه مطالعات محدودی درخصوص برآورد میزان نفوذ 
از بدنه رودخانه در حين عبور جریان سیل انجام شده است. همچنین 
این بررسی‌ها میزان نفوذ را تنها در روندیابی مستقیم جریان تخمین 
زده‌اند. لذا در رودخانه‌های فصلی و خشک که معمولا فاقد ایستگاه 
هیدرومتری می‌باشند. ضروری است جهت بررسی چگونگی طنیان یا 
فروکش کردن سیلاب میزان نفوذ از بدنه رودخانه نیز در فرایند 
روندیابی معکوس وارد شود تا بتوان با انجام محاسبات هیدرولیکی 
دقیق‌تر تراز آب یا دبی جریان را در نقاط موردنیاز به‌دست آورد. 
میزان نفوذ یا نشت به بستر رودخانه یا مسیل» ترکیب می‌شوند. سپس 
با روندیابی سیلاب در طول بازه و بررسی هیدروگراف‌های حاصل در 
نقاط مختلف» میزان نشت مشخص می‌شود (۸). 

بهرامی و همکاران (۳) به مدل‌سازی عددی نحوه پیشروی آب در 
جویچه‌های زراعی با روش روندیابی جریان و با مدل ماسکینگام- 
کونژ پرداختند که در آن معادله نفوذ لوئیس- کوستیاکف جهت برآورد 


خسروی و همکاران. توسعه روش روندیابی هیدرولیکی معکوس سیل در خشک‌رودها با درنظر گرفتن میزان نفوذ يا نشت ‏ ۸۰۷ 


میزان نفوذ استفاده شده بود. آن‌ها با مقایسه نتایج مدل عددی خود با 
خروجی نتایج نرم افزار 5180000 که روندیایی سیلاب را به 
روش‌های سینماتیک» اینرسی صفر و هیدرودینامیک مدل می‌کند» 
بیان کردند دقت روش ماسکینگام- کونژ به‌ویژه در جویچه‌های کوتاه 
و دبی ورودی بزرگ‌تر» بیش‌تر است. 

قبادیان و خلج (۱۶) مدلی عددی ارائه کردند که در آن معادلات 
سنت‌حونانت. با روش تفاضل‌های محنود حنل شند. رابطله تشست 
موسکات به گونه‌ای در معادلات مذکور گنجانده شده تابه طور 
همزمان نشت در مسیر جریان را نیز محاسبه کند. مدل ارائه شده با 
دقت مطلوبی هیدروگراف خروجی از بازه رودخانه زنگمار را پیش‌بینی 
نمود. محمدی و قبادیان (۲۱) مدل عددی را جهت تخمین مقدار 
نشت در رودخانه‌ها تهیه کردند که برای برآورد نشت رابطه گرین- 
آمپت در مدل مذکور با معادلات جربان غیرماندگار ترکیب شد. در 
مدل مذکور مقدار نفوذ در بازه‌ای از رودخانه قره‌سو در حدود سه 
درصد از حجم سیلاب ورودی به آن برآورد شد. 

نوردونین و همکاران (۲۴) تقریب موج پخشیدگی معادلات سنت- 
ونانت در روندیابی را با معادله نفوذ فیلیپ و جربان آب زیرزمینی 
ترکیب نمودند. خروجی مدل در مقایسه با داده‌های مشاهداتی» چشم- 
انداز امیدوارکننده را برای محاسبه مقدار نشت در جریان‌های ماندگار و 
غیرماندگار نشان داد. چنگ و همکاران )٩(‏ یک مدل روندیابی سیلاب 
شامل نفوذ از کف کال را پیشنهاد کردند. در این مدل» ییک روش 
شبیه‌سازی نفوذ از کف براساس تئوری نفوذ هورتون» توسعه داده شد. 
روش ماسکینگام-کونژ با جریان نفوذی از کف به عنوان خروجی 
جانبی» برای روندیابی سیلاب مورد استفاده قرار گرفت. نتایج 
شبیه‌سازی نشان داد که دقت این مدل بالاست و علاوه بر ساده 
بودن برای شبیه‌سازی و پیش‌بینی سیلاب در مناطق خشک و 
نیمه خشک مناسب است. 

محمدی و همکاران (۲۲) برای تخمین میزان نفوذ. مدلی عددی 
برای مقاطع نامنظم رودخانه‌ای تهیه کردند. در این مدل برای برآورد 
نفوف رابطه گرین-آمپت با معادله‌های جریان غیرماندگار چفت شد. 
ننایج این شبیه‌سازی نشان داد که وقتی از رابطه گرین-آمپت برای 
محاسبه نفوذ استفاده می‌شوده مقدار عمق نفوذ تحمعی و شدت نفوذ 
در هر یک از مقاطع در طول مسیر رودخانه قابل محاسبه و برآورد 
است. نویسندگان پس از مقایسه نتایج این مطالعه با تحقیقات قبلی 
بیان نمودند رابطه گرین- آمپت برای تخمین مقدار نفوذ یا نت 
نسبت به روابط دیگر از دقت بالاتری برخوردار است. پاتوواری و سارما 
(۲۵) یک مدل هیدرودینامیکی اصلاح شده جدید را با قابلیت تعیین 
دقیق مشخصات جریان در مناطق با نفوذپذیری بالا پیشنهاد کردند. 
این مدل با استفاده از معادلات جریان غیرماندگاره در ترکیب با معادله 
نفوذ گرین-آمپت به عنوان معادله حاکم» توسعه داده شد. برای حل 
معادلات حاکم» روش تفاضل محدود ضمنی به کار رفت. 


در مطالعه حاضر که به روندیابی معکوس هیدرولیکی سیلاب در 
گرین-آمپت در مدل عددی تهیه شده وارد شده است تا بتواند 
روندیابی معکوس را با دقت مناسب و در شرایط هیدرولیکی مسیل‌ها 
9 خشک‌رودها انحام دهد. بر اساس بررسی‌های انحام شدهء تاکنون 
هیچ پژوهشگری مقدار نفوذ را در روندیابی معکوس هیدرولیکی 
سیلاب مورد بررسی و تحلیل قرار نداده است و تنها توکلی و 
در حضور نشت يا نفوذ و به روش هیدرولوژیکی ماسکینگام غیرخطی 
به انجام رسانند. 


مواد و روش‌ها 
معادلات حاکم 

حرکت موج سیلاب در رودخانه‌ها از نوع جریان غیرماندگار متفیر 
تدریحی می‌باشد. در جریان‌های متغیر تدریجی شتاب قائم ناچیز و 
قابل صرفنظر کردن بوده» ولی اثر اصطکاک کانال قابل ملاحظه است 
و می‌بایست در محاسبات مدنظر قرار گیرد. معادله‌های حاکم بر 
جریان‌های غیرماندگار متغیر تدریجی به معادله‌های سنت-ونانت 


بوده و به‌صورت روابط زیر خلاصه می‌شوند (۱): 


(۱) پیوستگی: ۳ ۸ 90 
< ر 0 - ست + سک 
0 
2 
۳ ۳ 1 , 20 1 
۰۸1۸ * ۸ 


(۲) اندازه حرکت: 
(0< ( ,و و 8 

که در آن‌ها ۵ دبی جریان» ۸ سطح مقطع عرضی آبراهه 4 
جربان جانبی ورودی 8 شتاب گرانش» 5 شیب بستر و ٩۶‏ شیب خط 
انرژی ناشی از مقاومت جریان است. 

در تمام مدل‌های روندیابی هیدرولیکی جریان از شکل کامل 
معادله پیوستگی استفاده می‌شود. اما معادله اندازه حرکت در 
شکل‌های مختلفی که از حذف برخی اجزای آن به‌دست می‌آیند» مورد 
استفاده قرار می‌گیرد. ساده‌ترین مدل روندیابی هیدرولیکی, مدل موج 
سینماتیک است که در آن از اجزای مربوط به شتاب موضعی» شتاب 
انتقالی و نیروی فشاری در معادله اندازه حرکت صرف نظر می‌شود. 
یعنی فرض می‌شود ,5 وگ که بدین معنا است که نیروهای 
اصطکاکی و گرانشی با هم در تعادل هستند. در مدل موج پخشیدگی 
از اجزای مربوط به شتاب موضعی و انتقالی صرف نظر شده اما نیروی 
فشاری در آن لحاظ می‌شود. 


۸ نشریه آب و خاک. جلد ۰۳۵ شماره ۶ بهمن - اسفند ۱۴۰۰ 


1 1-۳۷۳۳ 


شکل ۱- شمای چهار نقطه‌ای با 7-0.5 (روش پرایسمن) 
(۱۱۵)۵۵0 یکت (۳) 05ص )۷1 عصصعمدامه )صتمه- ت۲۵ -1 مساعز۲ 


و در نهایت در مدل موج دینامیکی معادله اندازه حرکت به‌طور 
کامل و بدون ساده‌سازی در نظر گرفته می‌شود. به علت وجود اجزای 
خی فرسا ولا بسانم راما میتی سرام ای 
معادلات» به جز در مواردی که معادلات ساده شده‌انده وجود ندارد. لذا 
برای حل این معادلات باید از روش‌های عددی استفاده نمود. در این 
تحقیق از روش تفاضلات محدود. که روش غالب برای حل عددی 
معادلات جریان غیر ماندگار یک بعدی است. برای حل معادلات 
غیرخطی سنت-"ونانت استفاده خواهد شد. پرکاربردترین و قوی‌ترین 
روش به تایید عموم. شمای تفاضلی ضمنی چهار نقطه‌ای یا همان 
شمای پرایسمن است. از دلایل کاربرد گسترده این روش می‌توان به 
محاسبه همزمان دبی و عمق جریان در هر گره محاسباتی» امکان 
انتخاب گام زمانی بزرگتر و ضمنی و پایدار بودن روش اشاره نمود 
(۳۰. 

همان‌طور که در شکل ۱ نشان داده شده است. دامنه حل در 
شمای پرایسمن از شبکه‌ای به ابعاد (2,...,۷,/-) ۸:7 تشکیل 
شده است. تقریب مشتق زمانی و مکانی در نقطه م که در وسط بازه 
طولی :۸ است انجام می‌گیرد. این نقطه تنها می‌تواند در امتداد 
محور / یا زمان توسط تغییر مقدار ضریب وزن‌دهنده 0 جابجا گردد. 

اگر / یک متغیر وابسته مانند دبی جریان یا عمق جربان باشد. 
مقدار متغیر و مشتق‌های زمانی و مکانی آن بر اساس شمای چهار 
نقطه‌ای یا جعبه‌ای پرایسمن از روابط زیر محاسبه می‌شود: 


۳ زد م_ 1+[ اف 
) ۱ ٍ نکم ب ۵ 1+ (-1) نم 9 
(1> م > 0) ۸ ۸ 0 
۴ ۳ | 2 
ز دنلب ؛ 1+ ز (1-0) یمد 
(00۱) . ۸۰ ۸ 2 


(م -1)-+[ ۵(۲ ۵+0 |مع ,1 
(۵) [ بر( )+ ۵ ] 
فرد (۱۵) نشان داد که در حالت 7.0 > 0 > 0.3 شمای مذکور 
برای هر بازه زمانی به‌طور غیرمشروط پایدار است و مقادیر ۰/۵۵ تا 
۶ را برای آن پیشنهاد کرد. با جای‌گذاری ترم‌های تفاضل محدود 
حاصل از روابط بالا در روابط ۱ و ۲ و بعد از ساده‌سازی معادله 


پیوستگی و اندازه حرکت گسسته‌سازی می‌شوند. 


روندیابی معکوس به روش موج سینماتیکی 

همانطور که ذکر شد در مدل موج سینماتیکی گرادیان فشار و 
شتاب در مقایسه با ترم شیب کف کانال ناچیز است و نیروهای 
اصطکاکی و ثقلی با یکدیگر در تعادل‌انده پس حرکت موج را بیشتر با 
معادله پیوستگی توضیح می‌دهند (۱۰). 
می‌پردازد و به علت عدم نیاز به شرایط مرزی پایین‌دست بسیار 
پرکاربرد است. از آن‌جا که در اين روش معادله پیوستگی به‌طور کامل 
استفاده می‌شود. لذا معادله اندازه حرکت یعنی رابطه ۲ به‌صورت رابطه 
زیر خلاصه می‌شود: 
(۶) ,2 وگ 


رابطه فوق را برای استفاده در معادلات تفاضلات محدود می‌توان 


به فرم زير بازآرایی نمود: 
)۳ ۲*۶ 2 ۸ 


که ضرایب 0 و / با توجه به معادله جریان مورد استفاده قابل 
مقدار 5 و 0 برابر خواهد بود باء 


خسروی و همکاران. توسعه روش روندیابی هیدرولیکی معکوس سیل در خشک‌رودها با درنظر گرفتن میزان نفوذ یا نشت ۰ ۸۰٩۹‏ 


)۸ 25 
۳3 ۳ 
متغیر ۸ یا سطح مقطع جریان را می‌توان با مشتق‌گیری از رابطه 
۷ حذف کرد: 
)8 ۵و 9۸ 
9 2 "0 


و رابطه بالا را در رابطه پیوستگی جای‌گزین نموده تا معادله موج 
سینماتیک ناشی از تغیبرات 0 به‌دست آید: 


20 99 ۱۰ 
2 + 0/0 ۴ ۳ 

معادله بالا را می‌توان خلاصه‌تر نیز نوشت: 

80 2۵ ۳ 
0۲ 

کد: 

1 11 

(۱۲) اج 
۳0۵ 


سرعت موج سینماتیک است. با جایگذاری مشتقات مکانی و 
تفا ال شاه رایس ی را نا ۵ فوساداتم ۱۱ 


خواهیم داشت: 
۱۳ 1 ۱ - 0 (0 0( ۱ 0 ۳ 
0 ۸ ۸ 
رب 9:1 0 
+ +[ 1 ۱ 0 3 
۱ ۳ 9 ۷ 


20 ,+ 
و رابطه بالا را بر حسب مقدار مجهول در حل معکوس یعنی 07 به- 
دست می‌آوریم: 


(۱۴) مس ما ور ی ارم ای و 
ور 2 ۳ ور 9 


11 


که ضرایب موجود در آن از روابط: 


۱۵ لش 
0 ۳ 


0۵۲ ۱۷ 


0۸ [ ف رن 


تشر ۸ 
به‌دست می‌آیند. در استخراج معادله ۱۴ علامت عبارت .4 نسبت 
به معادله پیوستگی اصلی تغییر یافته است. زیرا عبارت .4 در مدل 
روندیابی خشک‌رودهاء نشت يا نفوذ به بستر رودخانه را نشان می‌دهد 
و باید آن را با علامت منفی در معادله پیوستگی وارد نمود. 


روندیابی معکوس به روش موح دینامیکی 

در روندیابی موج دینامیکی هیچ گونه ساده‌سازی در معادله اندازه 
حرکت صورت نمی‌گیرد و معادلات سنت‌سونانت به‌طور کامل حل 
من ورگ 

در روندیابی مستقیم» با گسسته‌سازی معادلات پیوستگی و اندازه 
کت زا انشفاده از کر یب‌هاش م زنط ی مق قای سکتای وا 
شمای پرایسمن و اعمال معادلات در جهت محور مکانی برای هر 
یک از ۸7 المان مستطیلی (مطابق با شکل ۲ بین محور مکانی 
[< و ۷ تعداد ( 22۷-2 معادله به‌وجود می‌آید درحالی که تعداد 
محهولات در هر گام زمانی 27۷7 است. دو معادله اضافی مورد نیاز 
برای تکمیل معادلات موجود برای حل مسئله» توسط دو شرط 
مرزی یکی در انتهای بالادست و یکی در انتهای پایین دست بازه 
رودخانه یعنی در ۲-0 و 2-۷ تامین می‌شوند. شرط مرزی بالادست 
را می‌توان هیدروگراف ورودی یا تراز سطح آب اندازه‌گیری شده و 
شرط مرزی پایین دست را شرایط جریان نرمال. جربان بحرانی و یا 
رابطه دبی اشل تعیین نمود. همان‌طور که از شکل ۲ مشخص است 
محاسبات از زمان اولیه یعنی 120 شروع شده و در جهت افزایش گام 
زمانی تا زمان 7-7 ادامه می‌یابد. مجهولات در شمای چپار نقطه‌ای 
پراٍیسمن دبی و تراز سطح آب در گام زمانی مجهول یعنی 
ی سس 

در روندیابی معکوس که هدف تحقیق حاضر است با به کار بردن 
معلدلات پیوستگی و اندازه حرکت در جهت محور زمانی برای هر یک 
از 1-/۸ المان مستطیلی (مطایق با شکل ۳ بین محور زمانی 1-ز تا 
تعداد (2-/2/ معادله به‌وجود می‌آید ولی تعداد مجهولات در 
هر گام مکانی /21 است. دو معادله اضافی مورد نیاز برای تکمیل 
معادلات موجود برای حل مسئله. توسط شرایط مرزی در 0</ 
و 127 تامین می‌شوند: 


( ,۵( ) 9 
یا (۶) ره ( ۶ <) 7 


(۱۹) جریان یکنواخت قبل سیلاب 
(شرط مرزی در 120 , .> > 0) 


۰ نشریه آب و خاک جلد ۰۳۵ شماره ۶ بهمن - اسفند ۱۴۰۰ 


( ,۵ (۰) 0 
(۲۰) جریان یکنواخت بعد سیلاب یا( ,-(۸)۲۰ 


(شرط مرزی در 127 > > 0 


(۲۱) هیدروگراف و تراز آب مقطع (9۲-<(1) 0 


(شرط اولیه در اد ,17 > > 0) ( ۰۸۶( :)۱ 


در روندیابی معکوس» محاسبات از انتهای پایین دست یعنی .]2 
شروع شده و تا ابتدای بالادست مجرا یعنی تا 2-0 ادامه می‌یابد. 
مجهولات دراینجا ۳ ۱/۲۰۵۵ و ۳ هستند 

شایان دکر است گسسته سازی معادلات سنت-ونانت در هر گام 
مکانی منتهی به دستگاه معادلات 2۸ معادله و 2 مجهول می‌شود 
که از روش‌های تکراری مانند روش نیوتن- رافسون قابل حل است. 


60 , ۱ ۱ 


همان‌طور که دکر شد. در روندیابی هیدرولیکی سیلاب در 
خشک‌رودها که با لحاظ نمودن همزمان مقدار نفوذ یا نشت در بستر 
رودخانه است. معلدلات سنت-ونانت باید با ترکیب معادله‌ای که نشت 
را برآورد می‌کند به‌طور همزمان حل شوند. 

به منظور برآورد میزان نشت می‌توان از روابط متعددی که توسط 
محققان مختلف پیشنهاد شده است. استفاده نمود. این روابط جهت 
محاسبه دبی نفوذ به بستر کانال‌های خاکی یا محاسبه نرخ نفود از 
سطح زمین ارائه شده‌اند. 

به‌طور کلی مدل نفوذ انتخایی بهتر است علاوه بر درنظر گرفتن 
شرایط غیر اشباع ناحیه بین بستر و سطح آب زبرزمینی» تغیبرات غیر 
ماندگار عمق جریان بر روی میزان نفوذ را نیز لحاظ نماید شرایطی که 
تحقیق از روش گرین- آمپت برای محاسبه نشت استفاده شده است. 


۴ ۵0۸ ۱ 0 0۵ تب 
۸ ۱ 0۵( + ُ 


زر :1 + ۱ 


شکل ۲- شبکه محاسباتی اختلافات محدود ضمنی در روندیابی مستقیم سیلاب 
0۵۵ 1000 0۵6ص انامه تحصمتاهاه موی معصمع)ن ماتصت الم تامرصا -2 مساعز۳۲ 


۳۳9 5 مب 1 
۸( ۵۴ 0 ۵ 
۵۸۵ ۱ , () 2 ٍ 


شکل ۲- شبکه محاسباتی اختلافات محدود ضمنی در روندیابی معکوس سیلاب 
۵ 1000 موم۲ه۱ رز عارن ۳۵ افصمتاها میم معصمعع]0 ماتصظ اتمتامرصا 3 ماعز۲ 


در این روش نرخ نفوذ با گذشت زمان بر حسب پارامترهای 
هدایت هیدرولیکی اشباع» مکش و درصد رطوبت در جبهه مرطوب 
پش‌بینی می‌شود. ضمن آن‌که ارتفاع آب جمع شده بر روی سطح 
زمین نیز در محاسبه نفوذ تاثیر می‌گذارد. روش گرین-آمپت بر مبنای 
توسعه یک تئوری فیزیکی تقریبی است که راه حل دقیقی دارد. گرین 
و آمپت تصویر ساده شده نفوذ را به صورت شکل ۴ پیشنهاد کردند. در 
این شکل جبهه مرطوب یک مرز تیز است که خاک با مقدار رطوبت 
0 در پایین مرز را از خاک اشباع با مقدار رطوبت در بالای مرز جدا 
می‌کند. پس از زمان / از آغاز نفوذ جبهه مرطوب به عمق ‏ نفوذ 
می‌کند و آب تا عمق کم ۸0 روی سطح خاک جمع می‌شود. معادلات 
نفوذ در این روش به‌صورت زیر است: 
۳ ۶ 
20۷ 


۳۳ 29۷ ۳۳ 
۳۲۲ 1 +1]1-- ز 


۲ )( < 1 + ۵۵۷ ۱011+ 


که 0-[۵0<1, 6 درصد رطوبت اولیه. [1 تخلخل خاک, ۶ 
سرعت نفوذ (سانتیمتر بر ساعت)» 7 نفوذ تجمعی (سانتی‌متر)» / زمان 
(ساعت)» ۸ هدایت هیدرولیکی غیراشباع (سانتی‌متر بر ساعت) و ۷ 
مکش رطوبتی خاک (سانتی‌متر) هستند. علامت منفی به این دلیل 


جایگزین ۷ در معادلات می‌شود. 


مدل عددی روندیابی معکوس 
۳ جهت روندیابی هیدرولیکی معکوس سیلاب به روش موج 
سینماتیکی و دینامیکی تهیه شد. در این مدل‌های عددی با ترکیسب 
رابطه نفوذ گرین- آمپت با معادله پیوستگی. محهولات عمق و دبی 
جریان در مکان‌ها و زمان‌های مختلف محاسبه می‌شود. همچنین دو 
مدل مذکور قابلیت تغیبرپارامترهای عددی و فیزیکی (پارامترهای 
هندسی و هیدرولیکی رودخانه و پارامترهای معادلات نفوذ) را در خود 
دارند و با دریافت این پارامترها به‌همراه هیدروگراف پایین‌دست جریان 
به‌عنوان ورودی» هیدروگراف جریان (مقادیر عمق و دبی جریان) را در 
انتهای بالادست رودخانه تولید می‌کنند. روند نمای نحوه انحام 
محاسبات در مدل معکوس روندیابی دینامیکی سیلاب در شکل 
شماره (۵) نشان داده شده است. 

جهت صحت سنحی نتایج مدل عددی توسعه یافته و بررسی 
قابلیت کاربرد آن. از مجموعه داده‌های اندازه‌گیری شده یکی 
هیدروگراف معروف لین (۱۰) و دیگری هیدروگراف حاصل از مدل 
آزمایشگاهی-صحراتی بمبشی جی (۳/ استفاده شده است که در 
شکل‌های ۶ تا ۸ مقطعم تیپ مجرا با رودخانه به‌همراه 
هیدرو گراف‌های مشاهداتی نشان داده شده است. سری زمانی داده- 
های اندازه گیری شده در مقاطع بالادست و پایین‌دست و همچنین 
داده‌های ورودی به مدل عددی شامل خصوصیات هنداسی و 
هیدرولیکی و خصوصیات خاک بستر نیزه به‌ترتیب در جداول ۱ و ۲ 
آورده شده است. 


۲ 
۱ عمق ب سطجنی 
۶ ۳۵۱۱۵ سطح‌زمین 
۳ هی 6:0۲ 
تاحیه مرطوب 
00 ۱۳4 
,1 
جبیهه رطویش 
۱ ۷۸ ۷۲۰۱۸2۵ 
سس سح 
4 4 ۱ ۰ 
۳ 9 
۰ سب | 
1 
ت ت 4 
لنه5 1(۳۷ ۶ 2 


شکل ۴- متغیرهای موثر در مدل نفوذ گرین - آمپت 


[06مح حمتاهاتصا امرجه هم صا مماطاه۲ ۲۳۲۵۵۵۷۲۵ -4 ۲12۱۲۵ 


۲ شریه آب و خاک. جلد ۳۵. شماره ۶ بهمن - اسفند ۱۴۰۰ 


حل ماتریس ضرایب 
۵ کایه101 0611 ٩۵1۷‏ 


۲ 4 و ۱ ؛ 
محاسبه دبی و تراز سطح آب 


بر آی مقطع (1-1-]) 
۱۲۵۱۰۲۹1۹۵ عابتا 
(1--[) ۲۵۸۵ 11۵۷۷ 20 [ع۲ع۱ 


محاسبه دبی و تراز سطح آب 


جدید برای مقطع (12-1) 
۵۵۵ ۱۷۵۱۲ قاتا 
(1-!(-[) ۲۵16 1100 0صه ۱2۲۵1 


۳۹ زمانی آخر 0 
0 ]1.25 
51003 


ورود بارآمترهای قیزیکی, عددی و 
شرایط اولیه و عرزی به مدل 
آ ۱ 2.0 اقتوزدام قعع۴ 
40 1۳10121 - ۲22۲2۲۵16۲5 
0۵۵ 000۳0۲۷ 


محاسیه پارامترهای هندسی و هیدرولیکی جریان در 
مقطع بابین‌دست ([(2[) 
عااسح ۳۲ 0صده عو0عع عاهاه‌اه) 


ع ۵۱ ۱0 عط) آ۵ 22۲27616۴5 
(1(حز) 690 واه 


برآورد مقدار نفوذ با رابطه گرین-آمپت 
عک ۲۵۸6 1۱۱۲۲۵۲۲۵۴۵ ۲۳5۱۲۱۵ 
/۸:۲۱۲/-۶۲۵۵۲) 1 


محآسیه ضرایب ماتریس محاسبهة ضرآیب ماتریس 


مجهول برای رابطه مجهول برای رابطه مومنتم 
و عاقابهآ) 
۲ 6 ۵0]111 12111۲ 
0 ۱۳۱۶ 
«متنعنوح 


۵ 24ات او ) 
۲ 002/1112۳016 022471 
نات ۷اناصتاصمت 1۳۱6 


شکل ۵- روند نمای مدل معکوس روندیابی دینامیکی سیلاب 
[۵0۵ عصتاه0: 1000 متصصهصرول م۲6۲۵ 0۴ افص ۲۱۵۷۲ -5 عاع۲۱ 


حل معکوس معادلات سنت-ونانت بد وضع ۲ است. در روندیابی 
معکوس جهت حل برخلاف حالت روندیابی مستقیم در راستای 
تقوبت موج پیش می‌رود و به این علت هرگونه خطا در مقادیر تراز 
سطح آب يا دبی جریان در طول محاسبات تشدید شده و منجر به 
ناپایداری در حل می‌شود (۱۲). اين ناپایداری عددی با ورود ترم نفوذ 
به معادله پیوستگی شدیدتر هم می‌شود. 

لذا در ابتدا نیاز است بازه مناسب پارامترهای وزنی شمای 
پرایسمن و گام‌های زمانی و مکانی طوری انتخاب شوند که مدل 
عددی علاوه بر همگرا بودن» هیدروگراف بالادست جریان را با دقت 
مطلوبی به‌دست آورد. بدین منظور بعد از تهیه مدل عددی رون دیابی 


1- 111-0 


معکوس و با استفاده از داده‌ها و يا هیدروگراف‌های مشاهداتی موجود 
آزمایش‌های عددی متعددی صورت پذیرفت تا گام زمانی و مکانی و 
محدوده مناسب پارامترهای وزنی ورودی به مدل جهت صحت سنجی 
آن استخراج شود. 

برای تحلیل حساسیت مدل نسبت به پارامترهای ورودی, با انجام 
آزمایش‌های عددی متعدد. بازه‌های پایداری مدل عددی در روندیابی 
معکوس سینماتیکی و دینامیکی بازای تغییر پارامترهای وزنی مدل 
یعنی 0 و م به‌دست آمد. در مواردی که به‌ازای ترکیب‌های مختلف 
پارامترهای وزنی» محاسبات قبل از رسیدن به مقطع بالادست واگرا 
شده و یا هیدروگراف خروجی بالادست دارای نوسانات مشخصی باشد 
آن ترکیب به عنوان بازه ناپایدار مدل روندیابی هیدرولیکی معکوس 
درنظر گرفته می‌شود. شکل‌های ۷ و ۸ نمونه‌ای از این ناپایداری‌ها را 
مه کت ات ات ای مها وت ان 
معکوس با شکست مواجه شده است. 


خسروی و همکاران. توسعه روش روندیابی هیدرولیکی معکوس سیل در خشک‌رودها با درنظر گرفتن میزان نفوذ یا نشت ‏ ۸۱۳ 


(الف) 


(ب) 


شکل ۶- مقطع تیپ رودخانه يا مجرا در زمان اندازه‌گیری, الف) رودخانه فصلی هاگز واش (داده‌های لین) و (ب) مجموعه آزمایشگاهی - میدانی 
(داده‌های بمبئی‌چی) 
زطءنه0عظ) [06مص اه)صمرتهمنی-6۱0 رط) ورح)ع0 عصصنا) ۱۲۵ جاوه ۲۷۲۷ عم‌طاهه۲۱ (ه) ممتامءی وومتن ۲۵ تن [مصصوول اهز۲۵] -6 ماهز۲ 


نوسانات ایجاد شده در نتایج آنقدر شدید است که در برخی از 
آزمون‌های عددی, برای دبی و تراز سطح آب مقادیر منفی حاصل و 
باعث واگرا شدن مدل روندیابی معکوس می‌شود. 

نتایج این قسمت از آزمایش‌های عددی نشان داد که مساله 
معکوس برای مقادیر مشابه پارامترهای وزنی که برای روندیابی 
مستقیم مناسب است. قابل حل نیست. یعنی بازای مقادیر > ۵0 >0.5 
0 و03 - . دلیل این تفاوت» چرخش ٩۰‏ درجه‌ای راستای حل 
معادلات سنت - ونانت در روندیابی معکوس می‌باشد (شکل ۲ و ۴ 
مشاهده گردد). در شمای پرایسمن بکار برده شده در روندیابی مستقیم 
پارامتر ۵ به‌عنوان پارامتر ضمنی‌ساز و باعث ایجاد اتلاف در محاسبات 
است. در مساله معکوس عملکرد پارامترهای وزنی برعکس بوده 
بنابراین پارامتر ۵ به‌عنوان ضمنی‌ساز شمای عددی تعریف می‌شود 
قح بایتدار شرا 0 عم و 0 خاضتا مس کرد سیون 
مشخص گردید هرچه مقدار م به یک نزدیک‌تر باشد. اتلاف شما 
فزایش می‌یابد که باعث حذف نوسانات ده و حلی پایدار را تضمین 
می‌نماید. در نهایت بازه‌های پایداری مدل هیدرولیکی معکوس برای 
داده‌های مشاهداتی موجوده در روندیابی سینماتیکی مقادیر > 0 > 0 
5 و 1.0 > > 05 و در روندیابی دینامیکی 0.6 > 6 > 03و 
0 > > 0.6 به‌دست آمد که با نتایج تحقیقات شیمکوویچ (۲۸ و 
۹ که بازه پایداری مدل روندیابی معکوس را 0.5 2 م و 0.5 > 0 
تعیین نمود تطابق مناسبی دارد. 

شایان ذکر است در روندیابی معکوس دینامیکی به‌علت پیچیدگی 
و غیرخطی بودن روابط. شدت نوسانات و ناپایداری‌ها نسبت به موج 


سینماتیک بیش‌تر است. در نتیجه. بازه پایداری پارامترهای وزنی مدل 
عددی در موج دینامیکی محدودتر است. جهت تعیین مقادیر مناسب 
گام زمانی و مکانی ورودی به مدل عددی و بالا رفتن دقت 
محاسبات» سعی بر آن شد عدد کورانت - نسبت سرعت انتشار حل 
تحلیلی به سرعت انتشار حل عددی- نزدیک به یک باشد. لذا ۸ و 
ال برای هیدروگراف بمبئی چی به‌ترتیب ۲ متر و ۵ انیه و برای 
هیدروگراف لین به‌ترتیب ۵۰ متر و ۲۰ ثانیه درنظر گرفته شد. 

مقادیر کوچک‌تر برای این دوپارامتر علاوه بر فزایش زمان انجام 
بالادست ندارد. 
طول مسیر قابل توجه بوده و نیاز است مقدار این تلشات در رونندیابی 
سیلاب اعمال شود. لذا برای نشان دادن اهمیت موضوع در ابتدا 
مدل‌های عددی توسعه یافته برای داده‌های موجود و در حالت نفوذ 
صفر ارو اجرا شد که خروجی‌های مدل در شکل‌های 9*٩‏ ۷۰ آورده 
شد. در این شکل‌ها هیدروگراف ورودی یا بالادست به‌دست آمده از 
روندیابی معکوس دینامیکی در هیچ‌یک از سری داده‌های در دسترس» 
با هیدروگراف‌های مشاهداتی تطابق مناسبی نداشته است؛ طوری‌که 
دبی اوج هیدروگراف بالادست در حدود ۲۵ تا ۳۰ درصد کمتر از حالت 


واقعی برآورد شده است. 


۴ نشریه آب و خاک. جلد ۳۵. شماره ۶ بهمن - اسفند ۱۴۰۰ 


جدول ۱- سری زمانی هیدروگراف‌های مشاهداتی 
کصامرهع0۳0 ود احممتاه ۲ عهودان ]۵ ومزدهء 11006 -1 12016 


هیدروگراف بمبئی چی 
جادرهع0۳0 ۱ تملآه حاصصوظ 
۰ ما ۰ 
دبی زمان دبی ِ 
۱ و ۵ 
خروجی خروجی ورودی ِِ 
۱۳۴0۳۲ نیت 1110۳ ۳ 
(1/5) 1206 (0ع8) 11۳006 ۰ (1/5) 12۸06 (نع) 
0 0 0 0 
12 15 13 4 
7 12 118 7 
139 15 116 11 
185 86 292 18 
211 93 217 25 
240 100 30 32 
3215 118 119 45 
323.8 136 134 54 
3238 151 131 86 
229 161 9 105 
217 205 253 143 
115 240 2317 13۹4 
95 3201 165 190 
55 3294 116 226 
23 559 1.5 283 
14 681 52 373 
05 196 23 530 
02 89 14 146 
0 953 0 916 


با محاسبه حجم زیر نمودار هیدروگراف‌های حاصل؛ حجم 
هیدرو گراف روندیابی شده بالادست در حدود ۴۰ درصد کمتر از حجم 
سیلاب واقعی به‌دست می‌آید که بیان‌گر آن است که برای بررسی 
طغیان يا فروکش کردن سیلاب باید میزان نفوذ از بدنه رودخانه نیز 
در فرآیند روندیابی معکوس وارد شود تا بتوان با انجام محاسبات 
هیدرولیکی دقیق‌تر تراز آب یا دبی جریان را در نقاط موردنیاز بهدست 
آورد. 

در مرحله بعد با وارد کردن ترم نفوذ از طریق رابطه گرین- آمپت 
در مدل و با انتخاب پارامترهای وزنی در محدوده‌های پایداه دقت 
مدل‌های عددی در روندیابی معکوس سینماتیکی و دینامیکی بررسی 
شد. از آنحا که مسیل‌ها قبل از وقوع سیلاب خشک بوده‌اند بنابراین 
درنظر گرفتن عمق جریان یا دبی برابر با صفر به‌عنوان شرط اولیه 
جلوگیری از کاهش دقت مدل و واگرا شدن آن» عمق جریان در تمام 
طول مجرا برابر یک سانتی‌متر به‌عنوان شرط اولیه درنظر گرفته شده 
است. 

هیدرو گراف‌های خروجی مدل عددی در مقطع بالادست در 
مقایسه با داده‌های مشاهداتی در شکل‌های ۱۱ و ۱۲ نشان داده شده 


هیدروگراف لین 
جارهع0۳0 ۲۱۲ 12۳6 
۰ ما ۰ ۰ ما ۰ 
دبی ت دبی ورودی 1 . 
خروجی ۱ ۰ 6 10۱/10۱۲ 
۲ 6 ۱0/10۷۲ 3 "0 
(د/(هه) 1266 0۵00۲ 500) (۳0/9) 56۵) ۳0 

0 0 0 0 

0 32092 0 996 
46 265 کل 1117 
83 2914 5.9 1190 
184 2360 117 1360 
189 4081 226 1652 
182 4300 284 1895 
15.9 4445 292 1919 
136 4640 3209 2089 
122 5 311 2102 
85 3247 271 2129 

2721 19 5490 17 
32012 13.8 53709 58 
32644 53 6291 34 
4057 11 6534 21 
178 21 6972 16 
43 29 71190 13 
5363 1.5 7360 017 

52303 0.6 9096 0 

15306 0 10۳1317 0 


است. به‌منظور بررسی میزان کارایی مدل عددی در تخمین 
هیدرو گراف بالادست از شاخص‌های میزان خطای نسبی و ضریب 
نش- ساتکلیف استفاده شد؛: 


100 ۳( 


1015 )( ۱ 1 


0 


۳ 2 


‌ ظ‌ 5 م0 2 


۷5 -1- +2 )۲۵( 


(0 10 01 ۱ 


7 <[ 


که مقدار دبی اوج» زمان اوج و حجم هیدروگراف بالادست 
محاسباتی. ,[مقادیر مشاهداتی این پارامترهاه ,رم « ور م6 به‌ترتیب 
دبی اوج و دبی متوسط هیدروگراف مشاهداتی و مر م6 دبی اوج در 
هدر رکفت معا یبا اس 

در مواردی که طراحی سازه‌های مرتبط با فعالیت‌های مهندسی 
رودخانه مدنظر است. مقادیر دبی حداکثر سیلاب آهمیت پیدا می‌کند. 
همچنین حجم زیر هیدروگراف که معرف حجم سیلاب می‌باشد در 
تعیین حجم مخازن ذخیره و همچنین محاسبه میزان تلفات جریان در 


خسروی و همکاران. توسعه روش روندیابی هیدرولیکی معکوس سیل در خشک‌رودها با درنظر گرفتن میزان نفوذ يا نشت ۰ ۸۱۵ 


بازه عبوری از رودخانه کاربرد دارد. 

در جدول ۳ و ۴ مقادیر دبی اوج. زمان دبی او حجم 
هیدروگراف‌ها و پارامترهای آماری حاصل از روابط ۲۴ و ۲۵ آورده 
شده است. هرچه درصد خطای نسبی کمتر و معیار نش- ساتکلیف به 
عدد یک نزدیک‌تر باشد نشان‌دهنده تطابق بهتر خروجی مدل عددی 
با داده‌های اندازه‌گیری شده می‌باشد. 

همان‌طور که از نتایج مشخص است هر دو مدل روندیابی 
معکوس هیدرولیکی استفاده شده با اعمال ترم نفوذ در معادلات با 


تفاس روگ اش رالاس زر رای موی ها سس 
داده‌های هی درو گراف بمبتی جی (هیدروگراف اول) دارای دقت 
بیشتری نسبت به هیدروگراف لین (هیدروگراف دوم) است. به‌عنوان 
مثال دبی اوج محاسبه شده به روش موج دینامیکی برای داده‌های 
هیدروگراف دوم دارای ۷ درصد خطا نسبت به مقدار مشاهده‌ای آن 
است در صورتی که برای داده‌های هیدروگراف اول خطای نسبی ۲ 
درصد به‌دست آمده است. 


جدول ۲- داده های استفاده شده در مدل عددی (۴ و ۲۰) 
(20 و4) ۱۵061 ایزهصصیه عمط سا 0عویا )۳2 -2 12016 


نوع هیدروگراف 
86 0۳0۵۲۵0۳ ز۲۱ 
(بمبتی چی) (لین) 
(نطمنهحاصوی   -‏ (مصصل) 
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شکل ۷- مثالی از ناپایداری مدل عددی روندیابی معکوس بازای 920.8 و 0.6-- هیدروگراف بمبئی چی 
احرع0۳0 و نطته‌حاصوظ -0,6عي و 00.8 10۲ واتازصاهاعصا 061مصه ممه۲ع۳ عظ) 0۶ عامرصهی م۸ -7 ۲12۱۲6 


۶ نشریه آب و خاک. جلد ۰۳۵ شماره ۶ بهمن - اسفند ۱۴۰۰ 
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شکل ۸- مثالی از ناپایداری مدل عددی روندیابی معکوس بازای 9-08 و 0-0.6- هیدرو گراف لین 
و00 وط مصصنا -0,6حم و ۵20.8 1۵۲ و)لازصاه)کص ۵0ص عوع ۲۵۲۷ 66 ۵۶ مامرصوی ص۸ -8 ۲12۷۲۵ 


0۱5 
020 ۱۵۱۱۵۱۷ قاتا ۰ 0.045 
)۵ 9۵۵۲۷() عه حه حه 00۳ 
۷ 1۳56۳۷0 4 کون . ٩۰‏ 
9 ض 5 
+ 6108 ۳ 3 
۰ 0,025 1 4 
+ 002 3 , 
2 2 
۰ 1 1411 گّ ۳ 
۳۹15 2 2 
1111 تك 


0۳۳0۱05 


۱000 800 0 . . 400 200 0 
زمان (ثانیه) 


11006 )۹۵۵( 


شکل -٩‏ خروجی مدل روندیابی معکوس دینامیکی بدون منظور کردن نفوذ - هیدروگراف بمبئی چی 
حارمع0۵ بط نمل‌نم‌ماروظ - ممناهماصا مص امه امد مه متام موه ۱۵۲ متصصحقطول ۵ موجن -9 مساوز۲] 
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شکل ۱۰- خروجی مدل روندیابی معکوس دینامیکی بدون منظور کردن نفوذ- هیدروگراف لین 
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خسروی و همکاران. توسعه روش روندیابی هیدرولیکی معکوس سیل در خشک‌رودها با درنظر گرفتن میزان نفوذ یا نشت ‏ ۸۱۷ 
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شکل ۱۱- مقایسه هیدروگراف بالادست روندیابی شده در مطالعه حاضر با هیدروگراف بمبئی چی 
ومد مصصنا از واه )همم مه 0مصنهاهان امرههم 0 «مالصا ما ۵ معتهم‌صم -11 و۲ 
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شکل ۱۲- مقایسه هیدروگراف بالادست روندیابی شده در مطالعه حاضر با هیدروگراف لین 
دروم0 و مصنا ۱ ای کمن عمط صز 0عصتهاهان حامرحهم بط ملک ما ۵ جمکزتوم‌صصمن -12 عساع۲1 


مقدار ضریب نش- ساتکلیف نیز برای هیدروگراف اول ۰/۸5 همان طور که در قبل ذکر شد و بر اساس شکل‌های ۱۱ و ۰۱۲ 
محاسبه گردید که بیان گر دقت بالاتر خروجی مدل عددی نسبت به دقت مدل عددی رون دیابی معکوس در شبیه‌سازی هیدروگراف 
هیدروگراف دوم با ضریب ۰/۸۲ می‌باشد. در مورد زمان اوج» بیشترین بالادست برای هیدروگراف اول (بمبتی چی) نسبت به هیدروگراف دوم 
مقدار خطا مربوط به رون دیابی معکوس سینماتیکی داده‌های بالاتر است که یکی از دلایل اصلی آن طول کوتاه بازه بين بالادست 


هیدرو گراف اول و به مقدار ۱۵ درصد است. و پایین دست جریان است (جدول ۲). طی آزمایش‌های عددی که در 


۸ شربه آب و خاک. جلد ۳۵ شماره ۶ بهمن - اسفند ۱۴۰۰ 


این مطالعه انحام شد» مشاهده گردید که با افزايش طول بازه» احتمال 
ایجاد ناپایداری و یا خطا در محاسبات بیش‌تر می‌شود. 

بعنوان مثال با افزايش طول بازه روندیابی معکوس سیلاب در 
داده‌های هیدروگراف لین» محاسبات در روندیابی سینماتیکی حداکثر 
تا طول ۱۵ کیلومتر و در روندیبی دینامیکی تا طول ۸ کیلومتر به 
سمت بالادست جریان دارای پایداری عددی بود که به علت 
ساده‌سازی صورت گرفته در معادلات روندیابی سینماتیک» 
هیدروگراف در طول بیش‌تری در جهت بالادست جریان روندیابی 
گردید. 

مقدار تلفات انتقال واقعی داده‌های در دسترس با مقایسه ححم 
زیر نمودار هیدروگراف ورودی و خروجی مشاهده‌ای قابل برآورد است. 
به‌ترتیب ۲,۱۲ و ۲۰۰۰۰ متر مکعب است که به‌ترتیب ۲۳ و ۲۴ درصد 
از حجم جریان ورودی به بازه‌های رودخانه یا کانل را شامل می‌شود. 

با بررسی خروجی نتایج مدل عددی و مطابق با جدول ۲ و ۴ 


دقت مدل‌های روندیابی معکوس در برآورد حجم هیدروگراف 
بالادست و در نتیحه میزان تلفات انتقال به‌ویزه برای داده‌های 
هیدروگراف اول مطلوب بوده و بیش‌ترین میزان خطای نسبی در 
برآورد تلفات باتوجه به نتایج به میزان ۲۵ درصد در مدل روندیابی 
معکوس سینماتیکی و داده‌های مربوط به هیدروگراف دوم به‌دست 
آمده است. در بین دو روش روندیابی معکوس هیدرولیکی, نتایج روش 
موج دینامیکی در اکثر آزمایش‌های عددی دارای دقت بالاتری بود» به 
این دلیل که در این روش معادلات سنت-ونانت به‌طور کامل حل 
مجرای مربوط به آزمایش‌های بمبئی چی, اثرات حذف ترم‌های 
موجود در رابطه مومنتم جریان در فرایند روندیابی سینماتیک 
محسوس نبوده و باتوجه به ارقام جدول ۴ ننایج هر دو روش 
روندیابی معکوس تقریبا نزدیک به هم می‌باشد. 


جدول ۲- دقت مدل‌های عددی روندیابی معکوس در برآورد پارامترهای شساخص در هیدروگراف بمبنی چی 
جارعم و تزهحاصصوظ ‏ ممامصصههمر ملامصتصصماع0 مصتهای ص کام0 مه آهمزم مه عصتت0 عو۵ ۲۵۲ 0 م۸ -3 12016 
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مدل معکوس دینامیکی مشاهداتی 
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جدول ۴- دقت مدل‌های عددی روندیابی معکوس در برآورد پارامترهای شاخص در هیدروگراف لین 
جوم و مصقا مز ممام‌صهتوم ۲ امصنصهاع0 مصاههتای 1 دمص امه ممابم عفه ۲۵۲۷ 0۶ رد۸ -4 12016 
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به‌طور کلی در مدل موج سینماتیک در صورت حذف اثر مقدار 
نفوذ بر روی حجم هیدروگراف فرونشینی موج اتفاق نمی‌افتد و 
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هیدروگراف پایین دست بدون تغییر شکل به مقطع بالادست منتقل 
می‌شود. از آن‌جا که فرونشینی امواج سیلاب در هر آبراهه به خوبی 


خسروی و همکاران. توسعه روش روندیابی هیدرولیکی معکوس سیل در خشک‌رودها با درنظر گرفتن میزان نفوذ یا نشت ۰ ۸۱٩‏ 


مشهود است. در نتیجه تعیین درجه اهمیت شیب بستر نسبت به بقیه 
جملات معادله اندازه حرکت که در مدل موج سینماتیک حذف 
می‌شوند. ضروری است. با این تفاسیر می‌توان دریافت که در مجاری 
با شیب زیاد. اهمیت جمله م5 از بقیه جملات بیشتر است. لذا استفاده 
از روش موج سینماتیک در صورت زیاد بودن شیب بستر پیشنهاد 
می‌شود. البته با افزايش شیب کف و در نتیجه سرعت جریان» جریان 
به سمت حالت فوق بحرانی تغییر کرده و در این‌حالت به دلیل حرکت 
مواج به سمت بالادست جریان؛ روندیایی معکوس سیلاب قاببل 
کاربرد نخواهد بود. 


نتیجه‌گیری 

در پژوهش حاضر دو مدل عددی موج سینماتیکی و دینامیکی 
توسعه یافته است که با استفاده از حل معکوس سنت-ونانت و کوپل 
کردن معادله نفوذ گرین- آمیت با معادله پیوستگی جریان» رون دیابی 
معکوس سیلاب در خشک‌رودها را همزمان با برآورد نفوذ در مسیر 
رودخانه انجام می‌دهند. برای صحت‌سنجی مدل در شبیه سازی 
سیلاب خشک‌رودها از دو مجموعه داده‌های اندازه‌گیری‌شده به 
نام‌های هیدرو گراف لین و داده‌های آزمایشگاهی- صحرایی بمبثی 
چی (سال ۲۰۱۲) استفاده شد. نتایج این مطالعه نشان داد که 
مدل‌های عددی مذکور توانایی شبیه‌سازی معکوس جریان را در 
مدل روندیابی معکوس هیدرولوژیکی ارائه شده توسط توکلی و 
همکاران (۲۷) امکان برآورد میزان نشت با استفاده از روندیابی 


معکوس سیل بررسی شده است ولی روش هیدرولیکی استفاده شده 
در تحقیق حاضر به‌علت حل کامل معادلات سنت-ونانت دارای دقعت 


منابع 


بالاتری نسبت به روش‌های هیدرولوژیکی است. همچنین معادله 
نشت استفاده شده در مطالعه توکلی و همکاران (۲۷) رابطة غیرخطی 
معروف به سینگ است که مقدار نفوذ در آن تنها تابعی از دبی جریان 
می‌باشد درصورتی که میزان نفوذ در خشک‌رودها علاوه بر شرایط 
هیدرولیکی جربان» به پارامترهای فیزیکی خاک بستر نیز وابسته است 
که توس راسله گزین > آمیت: رما دی اف زورفا رتیه 
شده است. 

میزان دقت مدل عددی توسعه داده شده در این پژوهش در 
برآورد دبی اوج هیدروگراف بالادست بالغ بر ٩۰‏ درصد با حداقل مقدار 
ضریب نش- ساتکلیف ۰/۸۲ و در محاسبه زمان اوج دقت مدل بین 
۵ تا ٩۷‏ درصد است. همچنین ححم هیدروگراف بالادست که نشان 
دهنده عملکرد مدل در برآورد نفوذ در مسیر جریان است با دقت ۷۵ تا 
۸ درصد برآورد شده است. محمدی و همکاران (۲۱) با جفت کردن 
معادله گرین- آمیت و معادلات سنت - ونانت در روندیابی مستقیم 
دینامیکی» میزان تلفات نفوذ را برای داده‌های هیدروگراف لین در 
حدود ۳۰۰۰۰ مترمکعب برآورد نمودند که ۸۰۰۰ مترمکعب بیش‌تر از 
مقدار واقعی آن (۲۰۰۰۰ مترمکعب) محاسبه شده است ولی در مطالعه 
حاضر میزان این تلفات در روندیابی معکوس دینامیکی ۱۹۰۰۰ 
مترمکعب برآورد شده است که نشان‌دهنده کارایی مناسب مدل عددی 
روندیابی هیدرولیکی معکوس می‌باشد. به دلیل رویکردهای متفاوت 
در محاسبه نفوذ (استفاده از روابط نفوذ در حالت اشباع و غیر اشباع 
خاک بستر)» روش‌های مختلف محاسبه تلفات در مدل روندیابی 
هیدرولیکی معکوس, می‌توانند تفاوت‌هایی را از نظر قابلیت کاربرد در 
روندیابی سیلاب و اثرات آن بر روی هیدروگراف خروجی مدل ایجاد 
کنند. این موضوع هم باید مورد توجه کاربران روندیابی معکوس باشد 
و هم برای تحقیقات آتی در این باره مد نظر قرار گیرد. 
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597-600۰ :(1003 ۱۵56۵2۲0 

۵۲ «162 مخ .صاملناه وعماموع۲ ماه ۷۷ .عصتاا۲۵ متحصقصل )تمتامرحصا صا 5126 ماد مصصتا ۵۶ ۳1]6۵)5 .1973 ,1.۱ ۳۲۲۵۵۵ 
338-۰ :(902 800120108وظ قعمعبا۴۵50 

همم عون ب۲۱۷۵ تمحصعم22 صا فعووما عمه6۵ع۹ 207ع)فصتا ۵۶ عصتامعصتاو۴ .2012 ۱۷۲۰ رهظ 20 ریبک صفتههصامطت) 
موه طفتاعحظ اه دوع 12) ,206-7214 :(3)6 معهصلهرنا مه ممتاهعتت ۵۶ متام صقتصه؟ ,۵061 

۸) ۵061 ۷۷۵۷۵ عمتفتلن چم 0عو2ه عصتانا10 ۲1۵00 عوز1۵۷۵ .2019 ۳۰ 4نطوهتطع۱۳ تمحصطه مه رطه ۵۵ ری هل ۴۱۵ 
۳۵۱۹128 18) .1003-1012 :(51)4 ط6270وم م2 ۷۷ همه 9011 ۵۶ لقصتاهل مفتصفط .رفصت ره ما۲ م۵ «لتتاد عقهم 
9اععتاوطاه طوناعطظ 

.عوعرظ 4محصطفه ۵۶ اهتنا توسمع۲ .عمتایای۱۷0 امصصجیت موم .2015 .ل تصصمطهترمه هه .6 نصته‌ووم۲۱ 
۳۵۲5121 

تون 0م)۲ع 1 ما از متام 11000 عورمب مک .2012 ,هه منوت همه ریگ فلل‌تامار ریک ۱۷۲2۸21 ,.بانظه فلووبامکز 
217-۰ :(12)1 9016668 موه طاعقظ مه وته۲22 لمتیتد۲ .عصفطهه عصتاتا۲0 ٩001226‏ 

اعصصحطه صععتاه احعصعمدمه مه عم دمتطدممتاهاه؟ انامانام .1 :1971 .1,0 لتمصیگ 4ص رب۷ صتلوزنا ولا مهم[ 
22-0۰ :13 ۳۲۵۲۵۱۵2 گه امصیامل ,موه 

۸) ۲۱۷6۲۵ 60۳0۵۴61۵1 طا 109965 ممتففتحصفصهه گه عصتامصصتافی 0ص عصتادا10 ۳۱۵۵۵ :2013 بک۱ صمیلدفامطت) مه رو تمصهطم]۱۷ 
0 مصوواء۲] مه طانه همتهعمتعمط! همه عصمجصه‌هممه]۱۷ م۴۱ وم عمصعهتصمت تقممتاه۲ 2۳ ,ویر ممفعتقطم :تاه عقعع 
(رصهلوهظ م) .صهی رصحتطه ] بمه0۲0 وش 

۲ هام60 اصدصه ۷ مج 0ج امحصش-جهععن) مصتام‌تامی .2015 5.۸۲۰ عنام‌مدگمطومک مه ریک صمتلدطمدت ریگ تمصهطم]۱۷ 
فصح وععصمزمگ ممتاموتتا ۶ه اقصتهل صهتصهیا .هبتر صا مصتاناهع ۶1000 هصتتل . وعقوم1. صمتفعتصعصه ۵۶ عصتاقصتادم 
(امعقتاوهاه حوتاعص طز صهتوع۴ 12) ,143-153 :(39)01 عطلمع1ع۲۸ 

۰ ۷ ۲۱۱2۵۱ 220 رب 9۵۵۷ .0.۷ 0۷۵ظ00عصه1 ,۷ 12600۷ وت قالتک ررت) 100حعظ رم صقطونا و1 01۵060 ر.ظ صتعم]۷ 
۶ هنال .هرطاتطصواظ عبر حاموتینک تیه لهتعصعطامه مط) هصمله مععقطهع عمکنامه امد صااه 0ص عصتانا۲۵ ۳1۵0۵0 
.262-۰ :368 

۵ ۳۹1۴۵۵01۴8۵ .2014 .با ۲۵/۵۲5 220 .۲.0 ۵01 هر طمامصانها! ما مع111 ,با 10ممصفطی رساه صزنالت۵۵ 
5 216۲ ۷۷ .0202 1۵7۵ ماه صتمتو ‏ مه ماه ممهاه . عصتود فامفصصقطل . موه اهمصمطاوه ‏ ص1۵ تا 
1474-1۰ :(5002 ۱۵562۲600 

۵۷1۵ ۲۱۷۵۲ 2 1 م۱ پرلمعافصن عصتادام؟ م۶ 061مص متحصحطول 0 بوط 011160مصظ خر .2017 ,کلیهر مهو 24 رو هماج 
60-۰ :(65)1 ومتصهه)معصهن0۲ ۲۱ مه زع۳۱۷۲۵۱۵ ۵۶ اقصتتاه؟ .عصم2 )ممحصلعنم 

مصصقهطه مهمن صا ۶109 لهماعصا مصتاهاناعع: 1۵0۶ محصمطمه لمزمصصبنه )تمتام‌حصا .2011 ,۳۱۰۸۷۲ تباید مه .۷۲:1 مفحصهطو 
دهع مه ۸1۵ ۲۷۷۲۲۱۵-۰۱5۰ .عممععومن روماممطهی ]ماه ۷۷ تمصمتامممنم] 15۳ 

م۷0۱0 2 وصلود قزر اهتعصمدامه ظا عنام ۶۱۵00 ووععم۴ .2021 .۱ تلهعدظ مصه ,۷ تامعصه2 ر.ظ تامعلله12۷ 
.صف رلفططمما که وتجهبنونا تمسملع۲ ,2021 ۳۵۵۰ 15-16 رقعصمع‌لممن ماما مفتمفی؟ ور ,مممتووه 
۳۵۲5121 

۰ ,۱۷۲۵0۱0۵۷ ۵ لمصتتاه .عصمتامنامع تصمصو۷ مصنجی مطا >م1 صماهامتم عوععص1 گم ممتانتآمو ,1993 ۴ 52۷۳۱۵۷۱۵2 
105-120۰ 

۵۶۲ 2۲6۷ ۱۲۲۵۵۵۲ طا عصتاتا10 11000 ووععصا م6 قاممصظ 0منام‌حصزه فط گه ممتاه‌تانمو۸ .2008 نک 52۷۳۵۷162 
.121-۰ :405 روعع61۵68 ۵۶ زصملمعه مفتاوط روعتورطنمع0 ۵ متتاتاوصا مطا ۵۶ فممتامءتاهظ .(0صق۵1ظ) 

۱۱6۷ .زممصتتمگ .ععتناهعبوط امصصقطه مهم۵ طا عصتاع0 مه لهمزعمصصی .2010 ,1 ,52/۴۵/1۵۷۷162 


4 
نشریه آب و خاک (علوم و صنایع کشاورزی) 8 آزمک 0صرح م۲۱۷2 ۵۶ مهم 
جلد ۰۳۵ شماره ۰۶ بهمن - اسفند ۴۰۰ ص. ۸۰۵-۸۲۲ باس 805-2 0۰ ,2022 ۷]21-۰-.۲60 و6 ,۷0 و35 ,۲۷۵ 


۵ ۲۵۱۱۵۱۱۵۲۵ رز ۱۷۱6600 مصت منم ۳۱۵۵0 ۲۱۲0۵۵۱6 رهم۱۹۵۲ 0۴ ۵۲۷۵۱0۵۳۱۵0۵( 
6 ۱۸/۱۱۵۰۵۱0۴ عمجم 


و ۹ ۱۹ 
04-12-01 :16061۷60 
05-01-72 :066060عظ۸ 


ما 21 وعتافتتمامهت2ط0 11000 مط غبامجاج 966060 وا ممتاعصصمصد موی ر0م1 2 تمه روععم «صمحظ هر مصمتا00 1۱06۲۵ 
۲ 200 ,1150109186 110 ز۵ تعزز مط صم مهو متتاممه۳۷0۲۵0 مظ و معقطا عمط ممتاهعم1 مل0عهو 2 ۵۶ دمن 
0 0۵4و مته فممتاهانهماقه بعصتادام 11000 ۲۵۷۵۲9۵ 1۱ .11000 معط ۵۶ محصتا فص 2 ۵0تناعحمطظ 0 و۳72 16۷۵ 266ناو 
0 0۲۳0006۵۵0 15 1۳0۵۲۵ ]علض مج بع2۳۲ مص ۵۶ ممتامعو صهه‌اومنا فص مع ماد 0 ومد ممتامعو صهه‌تاوه 0۵0۳ معط 
85 ۴ مصصتان۷ 26 روته1۱۷ 1 10005 انا بطاموهومط بط هموم هم مه متام اتمه متلنیمه بط مار م۵0 
۵۵ 21160 ,عصتان۷ ۱000 مد «مت)منالع؟ وتط1. .وعزه 0ج ۵60 فص مفص ممئلهت الصا ما مال وعقوعتعع0 ب«القتامه۲ع 
صرح له صا چته۲۱۷ تقصمف2ع؟ طاً هصتانام۳ 11000 ۶۵۲ داعم 06۷610۵60 ,۲۳۵۲۵10۲۵ .وههعه 2710 طا کصمعنصوله و1 روعووم1 
فص ,۳06685 عصتانام معط مد .ععقوما ممتفعتصصمصهط ۵ متقطتاوه ماق۵۲00۲و2 مه 0۲0۵۷108 10تامطو فصمتعع 86101-2110 
۸ 8 1۳0 وعوو10 م6۵۵۵ 0۲ طمتاهتصا عم غصتامععه مه مملاقناوع هچ اه 0مصتصاصمی مه فممتفنوه عمته0۲ع 
601 2 فصامهتع۳۷0۲0 عمتلنافی: مط مصتمتصین مه لهماصز معط عصملجة 11000 فط متام ها رفظ لقصقع ز۵0 
6 0 71 ولوعل طمتطان تاو تصوومتم مطا مر ,0مصتصصمامل وا وععوما ممتعفتصعصهع ۵۶ امه فطا رماصتوم 
0۳۵6۵۵۵ مطا رها 0عامصصتاوع عبر مطا معماه وععوما صمتاهعص1 ما رقوعته یه طا و00م1 گم مادام متانیجت بط 
م2000 اه عصتاناه 11000 و۲۵۷۵ موم م قع ۵و فلقمصظ. لهع1ععصصنته صا 06۲60قصمي فه9 وتطفممتاهام۲ 
فقط «لتتاد مظ که 50 ر۵ع۱۵۵۷/۱۵0 نان 0۶ ]65 ما 10 ,و1۷۵ آهتممممطاوه ۵ فممتنلمم متابهت بط تملصت هه پرمهتتاهمه 
۵۰ 1000 02116 زرط عوم۲۵ ز فعقوما صملففتصصعصهط 2021۷760 

اصجصه ۷ -صنعگ مطا عه م«مص! میج وسما ماصهته7 اتمه ممتو۲هع قفممتامتاوه 1 :ولمطع]۷ 0ص2 واهزعه)۱۷2 
عط رواممظ عصتبام 100 متابیهتبظ طا .عصمتاهتاوع مصتصعصصمصر فط صرح داوم فطع متاممد طمتطه رفممتامتاوم 
فصم1 میامزج۷ مد نومه فا ممتافتاوع مصتجممومصهر مصا تاه رلععتا ود ممتقناوع متام مطا ۵۶ م1 ماعامصمم 
۵061 ۷/۵۷۵ عتامصصمص! مط ور اع0مصظ عصتانامع ملابیهت بط )ومامرصصته م1 .عاممممم‌مهمی مصصمی عصدهمصصمر رها 0مصتقامان0 
مطا 1۱ .تن عتج مماجناوه حصتاصمجطمم فص مر م۵۲ معناووعزم 220 و0۲6۵ 121)تمص1 م6 ۲۵12660 فاممهمم‌صرمن فص ت۱۷ 
ماص حملها مد 10۲۵۵ مساوومتن مطا )ناه رلعااتصون مه وعم101 لهتاتمظ1 م6 ۲۵۱2060 فاممجمم‌صصمی فطا باعممظ ۷۷۵۷۵ صمتفتل 
۲ )ما۵۵ 106۲60قصهي و1 صمتاقتاوع حصتتصعمصطمصظر فطع رامق ۳۷۵۷۵ متصصتقصرول فطل صد تالمح فص صنامعمه 
۵ 10۳ ۹010۳۷۵۲۵ 2013 ظضش] ]۱۷۲۸ مد موز ۱۷۵۲۵ فهع 0۳0 ماهت2و۵و ,لاو کصمومزم مط ما .صماه) ماه 
6 عصتامنمه ها رقاع0مصه لهعزمصصتنه وففط) ما .فلمطامصه ۳۷۵۷۵ متصصقصرول مضه متمصصمصتکا ها متام 1000 متلنیهیل بط 
0 وهع0۱۸ )111676 طا مفهء 10 مج 06۵ مه بصمتامنامه متام فطا ط ممتمناوع صمتامع‌اصا امصشصوعین 
۰ 3۵) 0۶ 011۳600100 صروم‌طومنا فطع مر 40عاقلتاهم له 27 وعصصت 

ملقتله7ه 0 20 وام0مص همه 0۵۷7۵۱0۵60 فط ۵۶ ولناوعر معط ملقلتل۷2 مه تمهت ما خصمتفعیه عز 20 و6لناعع1 
4 5۵4 ۷۵۵ طموهع۲۵0ظ تدمزهماصصروظ مه طامهتوم ۲ مصقنا قه ط« مق قاق0 0متنافقممطر ۵۶ )6و 2 ,2001162011107 فا 
ط مووع00 ۲ صهه‌اوود مطا منم واممقور عصتانام 1000 متانیهته بط وووبع طاوها عقطا 07۵0و فلناوع۲ 
مصضا عمط صقطا مقنهععج معمصه ۳۵۲۵ 0202 طمهته00/ظ تطمزمصاصوظ فطا ۵ فافع مط .مجنمعه ماه20010۵۲۲ 
ما راصق ۲۵۲ .عفممتامعه مصههتامممل مه عاونا م۵۵ [۷2مص1 ما ۵۶ طاعجع1 تماعمطو فلا ما فتال طمقته070 ۶ 
جح 20 ع026 حامرحتع۳۷0۲0 مصمنا مطا ۲۵۲ 0مطاممهر ۱۵۷۵ مرول مطا رها 0ماهانملهه طامرهتعم0 بط اعلصا کن موتفطهفلل ۵2 
۰ ۷۷۵۵ 1202 طمهتع۷۱0:0ظ تطامزهحاحصدظ ۶۵۲ 0مصنهاداه تمه مط) ماتطاه رمتتلم۷ ۷60عوهاه فص مه ۵0عهم‌جط6۵ 79۵ ۵0۴ 61۲0۲ 
موتوه00/ظ مصوع‌موود ۵۶ ما۷ مطا عصتاقصتاوع ما قاملمصه متام مفتعبصا ۵ رمحتتاهعه2 فط. رععمامتعط [ 
.عاه0 محتع ملظ تدامزمماصحظ مطه ۶۵۲ «للهتممووع ماطامتنوع0 ۷/28 وعووما ممتعفتصعصهه ۵۶ امه معط «لتصمننمععومم 
عصضا 6 10 ۲۵12]60 0202 عط) مه عمط متامحصصمصت! ۲۵۷۵۲9۵ عطا 10۲ 2590 ۷28 وعووما مطتامططتاوه طا تمه ]فمطعنط 
۱( 

۵۵ 0۵1000 ۷۷۵۷۵ 0۱۱۵۵6 مطا ۵۶ ممصمممتنوم مطا رول‌مطاممه عصتانام 11000 متلنیهت0 بط مورم۲۵۷ و مطا ۳۲9۵۵ 
فتطا ما ۵۷۵0و «امامام‌مومه مته مممتلمنوه عصتصه7هع معط مه رقاممصصتهمون ناه دمم فط 107 ملقتتاهعة 20۲۵ 
۶ محممتانلجم لحنطمه مطا ماصعومتوع۲ من بجامرحتع۷0۲0 مصقنا مط مد ۵0مصنامصم01 102018 18 1016868 وتط1. ,۵00 


۶ ۵۱۱626 رعطنهءصلم۴۸ ممتامصصدامع مه ممتامعتس1 ۵۶ اممصاتهمعن1 ررمووع]۱۵ظ )002٩و‏ 220 )مهن ۲6,۲2۰ -2 20 1 
۰ (ر20 ۵0 ] ۵0 لوه۲ ۱۲ رومما0ومک تیدا هه نوخ 

(10]۵611 06۵ 110912۷1 .1.2 :معط مه م مت طممدمتزمن -۶) 

۷۷۵۸۱6۲ ۵22۷1 جهو0۲2طک1 ,میامن طمتهموم۴ 0منآ۸۵۵ ۵۶ ۲۱۵۵۵ رمصلزممصزعطظ 2۷۱1 طر ,۲.12 -3 

1۵۲: 3 6 


۲ نشریه آب و خاک. جلد ۳۵. شماره ۶ بهمن - اسفند ۱۴۰۰ 


۰ ۵۵0۵۱۵۲۵1 هه 

۶ 206 110۷ ۵21 مطا عصتاهصتاوه ما 90۶۵ ۷28 و۵0۵۱ لهعزتمصصتاه 0۵۷۵10060 فطا ۵۶ زمحتتاهم2 م۲ خمتفمومن) 
6 015002726۰ فنطا 40 1۵۱2)60 مصصتا مطا عصتمصتاوع ما 976۵ مه 850 م۵0 مضه بامهتع0۵ت بو جصهع‌تومدا فطا 
مطا طا مناج الصا مطا عصتامصصتاوع ما معصقصمون م20 مطا ومیهعن4ص فتاه رامعم بط فم‌تاومده عطع ۵۶ متام 
ماه واع0مص لهمزمصصنه مطا عقطا «مطه عاناوع: عفقط1 .رمهتتاهع2 98۵ ۲۵ 75 ۷108 ۵06160 مفله ۷28 بطم 110۷ 
مه وعووم1 860۵826 0۲ ممفلص29 م۵۲ رفته 1۳۷ آهمصمداوه ما پمقتتاهعع فاماهاومععه اه عصتاتام 11000 ۲۵۷۵۲8۵ 
۲ 01116160669 معلقصظط هه و0مط)مصظر معقطا روعقوم1 ممت)اهاللصد عصته‌لنملهه ما وعطمهممم2 اصمتعنل ما مبانز اصوع/نصوله 
۰ ۵2اه ۵ آنامتان ط۵ه۲ع10ظ عطا 


۷ 9/60 رعصتاتا10 6۷6۲9۵ رواتاتطحافصا [06مصظ لهع1رمصیااظ رماه ممتاهص1 رصمتامنامع امصمشهه۵ع0۳ :1۵۲۱۷۵۲9 
فصمحلنهملقه وماو 


